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ABSTRACT 

We investigate the effects of the convection in the hyperaccretion disk around a steUar- 
mass black hole, which is considered to be the central engine of gamma-ray bursts 
(GRBs), with simple analytical calculations. If the convective energy transfer in the 
vertical direction becomes efficient compared with the inward advective energy trans- 
port, the hyperaccretion disk is expected to be hotter and the neutrino emission due 
to the electron-positron annihilation becomes the most efficient cooling process. We 
find that the sequence of the thermal equilibrium solutions for the convective hyper- 
accretion disk would have a viscously unstable branch, especially when the viscosity 
parameter is relatively small (a < 0.01). This means that the sporadical mass ac- 
cretion onto a black hole would occur in this disk. We propose that this process can 
be the origin of the highly variable lightcurves observed in the prompt emissions of 
GRBs. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic explo- 
sions in the universe, which release large amount of energy 
up to 10^^ ergs in only 0.01 - 10^ seconds in the form of 
gamma-ray emissions, and they show violent time variabil- 
ities in their emissions {St ~ 1msec). Although the central 
engine of GRBs is totally hidden from our view, accord- 
ing to their enormous power, it is often argued that rela- 
tivistic phenomena should be involved. The most popular 
model for the central engines of GRBs is a hyperaccreting 
black hole model (Narayan et al. 1992, 2001). In this model, 
the energy of relativistic jets that makes intense gamma- 
ray emissions is produced via the accretion of a massive 
disk (0.1 — IMq) onto a stellar-mass black hole. Such sys- 
tems are expected after several energetic phenomena such as 
mergers of double neutron star binaries (Eichler et al. 1989), 
neutron star/black hole binaries (Paczyriski 1991), white 
dwarf/black hole binaries (Fryer et al. 1999), black hole/He 
star binaries (Fryer & Woosley 1998), and failed supernovae 
(or coUapsars; Woosley 1993; Paczyriski 1998; MacFadyen & 
Woosley 1999). Their mass accretion rate is supposed to be 
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as high as 0.01 — IOMqs^^, and then the accretion flows be- 
come so optically thick that photons are almost trapped in 
the flows. Therefore radiative cooling is not efficient in this 
case. Alternatively, since the density and temperature can 
become very high (p > lO^g cm-^ T > 10^°K), the cooling 
via thermal neutrino emissions will become efficient. Such 
disks (or flows) are called "Neutrino-Dominated Accretion 
Flows" (NDAFs). This neutrino emissions may provide the 
energy deposition enough to making relativistic jet by neu- 
trino annihilation above the disk (Popham et al. 1999). The 
steady structure of NDAFs has been studied by many au- 
thors (Popham et al. 1999; Narayan et al. 2001; Kohri & Mi- 
neshige 2002 (hereafter KM02); Di Matteo et al. 2002; Kohri 
et al. 2005; Gu et al. 2006: Chen & Beloborodov 2007; Liu 
et al. 2007; Kawanaka & Mineshige 2007; see also Chapter 
10.6 of Kato et al. 2008). As for their time-dependent behav- 
ior, some authors performed hydrodynamical simulations in- 
cluding neutrino cooling (Ruffert & Janka 1999; Janka et al. 
1999; Proga et al. 2003; Lee et al. 2004; Rosswog 2005; Fuji- 
moto et al. 2006; Setiawan et al. 2006; Nagataki et al. 2007; 
Shibata et al. 2007; Metzger et al. 2008; Sekiguchi & Shi- 
bata 2010; Carballido & Lee 2010; Taylor et al. 2011), and 
some authors investigated analytically the disk instabilities 
which may occur in NDAFs and originate the fast variability 
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in the prompt emissions of GRBs (Janiuk et al. 2004, 2007; 
Masada et al. 2007). 

In this study, we investigate the structure and stabil- 
ity of NDAFs by taking into account the convective en- 
ergy transfer in the vertical direction, which has not been 
considered in detail in the previous studies. In the con- 
text of standard accretion disks which is radiation pressure- 
dominated, some authors have inferred the importance of 
vertical convective motions (Bisnovatyi-Kogan & Blinnikov 
1977; Shakura et al. 1978; Milsom ct al. 1994; Blaos & 
Socrates 2001; Sg,dowski et al. 2011) and some recent stud- 
ies have shown the importance of convective motions in the 
hyperaccrction flows (Milosavljovic ct al. 2010; Sekiguchi & 
Shibata 2010). Due to the vertical convection, the internal 
energy stored in the gas would not be advected inward but 
transported upward, which alters the equation of energy con- 
servation. We present the equilibrium solutions of convective 
NDAFs, discuss their stability, and propose the new scenario 
for the evolution of NDAFs that can account for the violent 
time variability in the prompt emission of GRBs. 



2 MODELS 

In the following discussion, we adopt equations describing 
the steady structure of a hyperaccrction disk based on New- 
tonian gravity. Angular velocity of gas particles can be ap- 
proximated as Keplerian: f2(r) — {GM/r^)^^^ , where M is 
the mass of the black hole. Expressions for mass conserva- 
tion, angular momentum conservation, hydrostatic balance 
and a-viscosity are given, respectively. 



and the neutrino cooling rate Q is composed of four terms: 
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where Vr, S, f, p and H denote the radial velocity, surface 
mass density, kinematic viscosity coefficient, pressure and 
scale height of the accretion disk, respectively; M and a arc 
the mass accretion rate and viscosity parameter. Here the 
pressure p is composed of three terms: p = Prad -I- Pgas -I- Pd, 
where Prad, Pgas, Pd are the radiation pressure, the gas 
pressure, the pressure of the degenerate particles, respec- 
tively (for their mathematical description see Section 2.1 of 
KM02). We set inner radius of the accretion disk, nn to be 
the radius of the innermost stable circular orbit around a 
Schwarzschild black hole, rin = 3rs, where rs ~ 2GM/(? 
is the Schwarzschild radius. In addition, we should consider 
the equation of energy conservation: 
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where and are the viscous heating rate and neutrino 
cooling rate per unit surface area. As for the right hand 
side, the energy dissipation rate per unit surface area can 
be described as 
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rems "r '/plasmon; 

where qNe, Qe+e- i ^brems and ^piasmon are the neutrino emis- 
sivity due to the electron/positron capture by a nucleon A'^, 
due to the electron-positron pair annihilation, due to the 
nucleon- nucleon brernsstrahlung, and the plasmon decay, re- 
spectively (see Section 2.3 of Kohri et al. 2005). The left 
hand side is often called as the advective cooling rate: 
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which means the inward flux of the gas energy along the flow. 
Here s denotes the entropy per unit mass, s = (srad-l-Sgas)/p, 
where Srad and Sgas are the entropy density of the radiation 
and the gas, respectively (see Section 2.3 of KM02). In the 
following discussions, we approximated ds/dr as s/r, and 
define the total cooling rate as = + Q^^^- It would 
be obvious that the radiative cooling is negligible. 

In order to see effects of convection in an hyperaccrction 
disk, we should estimate a typical timescale for the convec- 
tive motion of the blob in the vertical direction tconv and 
compare it with timescale for the inward advection tadv. 
The acceleration of the blob in the accretion disk due to 
the vertical buoyancy would be the same order of magni- 
tude as that due to the gravitational force, i.e. ~ Q{r)'^z 
where z is the vertical coordinate measured from the equa- 
torial plane. Then the convective speed along the vertical 
direction comming from Zmin (<SC z) should be 
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and then the timescale for the blob which is emerged at 
z{< H) and transported convectively can be estimated as 
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On the other hand, advection timescale can be esti- 
mated as 
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Here we concentrate on cases with M < 10~^MQsec~^. Ac- 
cording to KM02 etc., such an accretion disk is very optically 
thick with respect of photons and not so dense or hot enough 
to emit neutrinos efflcicntly. However, when the convection 
timescale is shorter than the advection timescale, the en- 
ergy dissipated in the accretion flow would not be advected 
radially inward, but would be convected upward, and as a 
result the gas at a certain radius would be heated more efB- 
ciontly than in the case without convection. Let us evaluate 
structures of such accretion flows with simple analytical dis- 
cussions. The convective motion can prevent the advective 
energy transport when fconv < *adv, which would be fulfilled 
when z > Zc where Zc should satisfy 
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Therefore, the equation of energy conservation should be 
modified by multiplying the factor Zc/H to the left-hand 
side of Eq. ([5]) to be 

TE^Vr^^Q+-Q-. (13) 

Other effects should be much milder than this suppression 
as will be discussed later. 

Because of this suppression of the advective cooling, 
the significant amount of internal energy of the gas stay at 
the same radial position in the accretion disk. As a result 
the disk temperature would be much higher than in the case 
when the convection is neglected, and then the disk structure 
is expected to be significantly changed. In order to see this, 
we plot the sequence of the thermal equilibrium solutions 
for convective NDAFs on the (E, M) and (E,T) plane for 
various values of E, fixing the radial position r. 
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3 STABILITY OF CONVECTIVE NDAFS 

By finding the sets of physical parameters which satisfy the 
fundamental equations for the disk structure, ([1]), ([2]), ([3]) 
and ([SI), we can derive the sequence of the thermal equi- 
librium solutions. In Fig. 1 we show them on the (E,M) 
plane for r = 4rs, varying the viscosity parameter a. Here 
the mass accretion rate is normalized by the critical mass 
accretion rate. 



Merit = lel/Edd/c , 

where LEdd is the Eddington luminosity, 

-t/Edd = 47rGA/BHmpC/(TT, 



(14) 



(15) 



with proton mass rrip and the Thomson cross section ot — 
6.6 X 10~^^cm^. We can see that for lower values of a the 
thermal equilibrium curves have a part with negative gradi- 
ent, i.e. 



V dE 
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which means that this accretion fiow is dynamically unsta- 
ble (see, e.g. Lightman & Eardley 1974; Kato et al. 2008). In 
such accretion disks, we can expect the sporadical mass ac- 
cretion onto the central black hole, which may occur in dwarf 
novae and some microquasars with a radiation pressure- 
dominated accretion disk. This can be the source of violent 
time variability in the GRB prompt emissions with the mass 
accretion rate of a few times 1O~*M0S~^, which corresponds 
to 10^" - lO'^^ erg s-^ 

We also show in Fig. 2 the thermal equilibrium curves 
on the (E, T) plane with the same parameter sets as in Fig. 
1. Here we can see that in the whole parameter space the 
thermal stability criterion. 
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fdQ- 



dT J ^ ^ \ dT 
is satisfied, and so this accretion flow is thermally stable. 
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Figure 1. Thermal equilibrium curves of convective hyperaccre- 
tion flows on the (S, M) plane at r = Arg (solid lines) with 
a = 10~^ {green line), 10~^ {red line), 10~^ (cyan line) and 10~* 
{blue line), plotted from top to bottom. The central black hole 
mass is set to be 3Mq. The thermal equilibrium curves without 
taking into account the convection are also shown (dashed lines). 
Right regions with respect to the bending part are dominated by 
degenerate-electron pressure in which we are not interested in the 
current study. 
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Figure 2. Thermal equilibrium curves on the (S,T) plane at 
r = Arg with the same parameter sets as in Fig. 1. There is a 
notable difi'erence between the current case and the convection- 
inefficient one where the T — S relation did not depend on a in 
the normal slim-disk or NDAF solutions (see Table 1 of KM02). 



4 DISCUSSION 



Here we discuss the physical origin of the dynamical instabil- 
ity of NDAFs with vertical convective motions presented in 
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the last section. To sec this let us express the physical prop- 
erties of the accretion disk in terms of (E, M) or (S, T). As 
for the former set of variables, it is obvious from Eq. (6) that 
the viscous heating rate Q+ is proportional to M. We can 
also express with the latter set of variables as follows. In 
the accretion flows with a mass accretion rate much larger 
than Merit, the pressure and entropy are dominated by ra^ 
diation, i.e. prad > Pgas,Pd and Srad » Sgas- In this case, 
the scale height H would be proportional to T*E~^f2~^. 
Therefore, the viscous heating rate Q"*" w 2apH0, should be 
proportional to ar**S"^fi"\ 

On the other hand, as for the neutrino cooling processes, 
the emissivity by the electron/positron capture process in 
electron-nondegeneracy limit can be approximated as 

= 9-2 X 10- ergs cm- (^)' (^qJo^) , 
while the emissivity by the electron-positron annihilation is 

9e+e- = 4.8 X 10- ergs cm"^ s"^ ( lo^) ' (^^^ 

in the nondegeneracy limit. In the unstable branch found 
in Fig. 1 (and the corresponding part in Fig. 2), wc can 
see that the pair annihilation process dominates the elec- 
tron/positron capture process. Therefore Qv is approxi- 
mately proportional to T^H oc T"E"^ 

By equating the viscous heating rate and the neutrino 
cooling rate, Q'^ and , we can easily see that the tem- 
perature T is proportional to a^^^ M^^^^r~^^^° and does not 
depend on E, which is also found in Fig. 2, because both Q"*" 
and Q,7 are proportional to E~^ when described in terms of 
(E, T). Considering the same part of the thermal equilibrium 
curves on the (E, M) plane, as M oc f2"^Q+ oc T*E"\ we 
can see that M oc a"/*M"^/^°r^^''^°E"\ which is seen in 
Fig. 1 as a negative gradient part in each thermal equilib- 
rium curve. 

We can interpret the physical reason of the viscously 
unstable branch in the coiivcctivc NDAFs as follows; due to 
the vertical convection in the accretion disk, the advective 
energy transport would be significantly suppressed, and thus 
the disk temperature would become higher than in the case 
without the convection. In such a hot disk, electrons would 
not be degenerate and hence the electron-positron annihila- 
tion process would become an efficient cooling process. It can 
dominate the emissivity due to the electron/positron cap- 
ture by rmcleons (oc T^'p), which was the dominant neutrino 
emission process in the NDAF models without the convec- 
tion (e.g. see KM02). Therefore the temperature-dependence 
of the cooling rate would be changed and make the un- 
stable branch appear in the (E, M) plane. We can see in 
Fig. 1 that with smaller viscous parameter a the unsta- 
ble branch appears with smaller mass accretion rate M. 
This is because the advective energy transport is more in- 
efficient with smaller q, which makes the disk temperature 
high enough for the electron-positron annihilation process 
to become dominant even with smaller mass accretion rate. 

In particular, the typical timescale for this instability is 
the order of the viscous timescale, 




and therefore when r/H < 0(1), which is often realized 
in the innermost region of NDAF, this timescale is short 
enough to account for the time variability observed in the 
prompt emissions of GRBs. 



5 CONCLUSION AND SUMMARY 

In this study, we show that if the convective energy transfer 
along the vertical direction, which is naturally expected in 
the radiation-pressure dominated accretion disk, the ther- 
mal equilibrium solution for a hyperaccretion disk would 
have a viscously unstable branch. When this branch is re- 
alized, highly time-dependent mass accretion onto a black 
hole would occur. Such instabilities have been studied in 
the contexts of dwarf novae (Meyer & Meyer-Hofmeister 
1981; Hoshi 1979; Smak 1982) and microquasars which have 
radiation pressure-dominated accretion disks in themselves 
(Lightman & Eardley 1974; Shibazaki & Hoshi 1975; Pringle 
1976; Shakura & Sunyaev 1976; Abramowicz et al. 1988), 
and the time variable features observed in those sources 
are theoretically interpreted with the model of those in- 
stabilities. In the context of GRBs, the process discussed 
in this study may lead to the intermittent energy release 
from a hyperaccretion disk, and then the highly inhomo- 
gcneous jet would be launched, which can be the origin of 
the short-term variability in the prompt emission of GRBs. 
Judging from isotropic luminosities of observed long GRBs 
(Liso ~ 10^^~^^erg s~^), the accretion rate adopted here 
M ~ 10"^ - 1O"*M0S"^ ~ 10^°"^^erg s'^/c^ is relatively 
small. However, recent observations show that there is a pop- 
ulation of low luminosity GRBs (I/iso ^ 10''®^''® erg s~^; e.g. 
Galama et al. 1998; Campana et al. 2006), and the variabil- 
ities of their prompt emissions may be explained with our 
model. Our model may be also appropriate to describing the 
late time activity of GRBs such as X-ray fiares in afterglows 
(e.g. Burrows et al. 2005). 

We should note that the absorption optical depth for 
neutrinos is smaller than unity, according to Eqs. (55)- (57) 
in Kohri et al. (2005). This means that we do not have to 
consider the neutrino trapping effect in the current situation 
(T < 10"K, E < lO^Sg cm-2). 

We should mention the effect of magnetohydrodyanmic 
turbulence on the convection in our accretion disks. A hy- 
peraccretion disk should have a turbulent structure due to 
the magnetorotational instability (MRl), and it would help 
the mass accretion as the turbulent viscosity (cf. Balbus and 
Hawley 1991). Here, using a parameter, we can describe the 
turbulent viscosity i^t as Ut ~ aCsH where Cs and H arc 
the speed of sound and the scale height of an accretion flow, 
respectively. In our study a parameter is assumed to be relar- 
tively small value, < 10~^ , which may be due to the suppres- 
sion of MRI in a hyperaccretion flow (Masada et al. 2007). In 
other words, the typical scale and the speed of MRI turbu- 
lence in our accretion flows are much smaller than the scale 
height and the speed of sound, respectively. So, in our cases, 
we can say that the convective motion would not be dis- 
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turbed by MRI turbulence because the typical scale length 
and speed of convection are ^ H and ~ Cs, respectively. 

The convective energy transport along the vertical di- 
rection can also affect the structure of the slim disk, in which 
the optical depth for photon r is so large that the photon 
trapping is efficient. If the timescale for the vertical convec- 
tion is shorter than that of radiative diffusion, the photons 
generated in the equatorial plane can be convected upward 
to the disk surface and emitted out efficiently. This means 
that the effective optical depth would get much smaller 
than that in the case without convection. This modification 
should appear in the radiative cooling term of the energy 
conservation equation approximately by multiplying a fac- 
tor of Zc/H , to the optical depth, i.e., like Tcfi ~ r x Zc/H . 
This modification in r makes the equilibrium curve for the 
standard disk shift to a larger E and a higher T regions such 
as E ~ 10* - lO'^g cm"^ and T ~ 10*if , respectively. In ad- 
dition to these effects in the multi-dimensional treatments, 
we should also consider a diffusion of photon along with the 
z-sxis, which also affects the photon-trapping region (Kohri 
et al. 2007). Therefore although the current NDAF region 
with E ~ 10" - lO^g cm-2 and T ~ W° K discussed in 
this letter should not be affected by these modifications, to 
know the whole structure of the disks, we would have to per- 
form much more precise numerical simulations with reliable 
treatments of the convection, the diffusion and the neutrino 
transfer. 
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